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INTERFERENCE REDUCTION BY STEP FUNCTION 

REMOVAL 

TECHNICAL FIELD 

This invention relates to reception of a signal. 

BACKGROXJND 

A received signal may include a desired signal from a 
desired source along with one or more undesired signals, such 
as, for example, a noise signal (e.g., additive noise such as 
white Gaussian noise) from a noise source, and/or an 
interfering signal from an interfering source (e.g,, main- lobe 
or side- lobe energy of the interfering signal) . The received 
signal also may include an offset component such as, for 
example, a DC (direct current) offset component, that may be 
undesirable. The offset component is an additional additive 
term and may be a constant offset such as, for example, a DC 
offset, or may be a non-constant offset such as, for example, 
a step function. 

To extract the desired signal from the received signal, 

characteristic parameters (e.g., data bits, frequency offset, 

DC offset) that model the received signal may be estimated. 

It may be desirable to perform preprocessing of the received 

signal prior to estimating the characteristic parameters, such 
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as, for example, estimating the offset and removing its 
effect. For example, the offset may be estimated as a mean of 
the received signal and the mean may then be subtracted from 
the received signal . 

The offset of the received signal may vary, for example, 
because of variation in the interfering signal. Such 
variations may cause the mean of the received signal to 
provide a poor estimate of the signal offset. Subtracting a 
poor estimate of the offset would then bias the estimates of 
the characteristic parameters and lead to inaccurate results. 

DESCRIPTION OP DRAWINGS 

Fig. 1 is a schematic diagram of a communication system 
configured to estimate and correct a signal having an offset 
that may be modeled as a step function. 

Fig. 2 is a diagram illustrating two time -division, 
multiple-access (TDMA) users that are not aligned in time and 
that may cause interference which may appear as an additional 
signal offset in the form of a step function. 

Fig. 3 is a schematic diagram of a receiver that may be 
used with the communication system of Fig. 1. 

Fig. 4 is a diagram illustrating a step function that may 
be used to model the signal offset of the system of Fig. 1. 



Fig. 5 is a schematic flow diagram illustrating a 
systematic process for offset correcting a signal having an 
offset that may be modeled as a step function by removing the 
step function of Fig. 4 from the received signal. 

Like reference symbols in the various drawings indicate 
like elements. 

DETAILED DESCRIPTION 

For illustrative purposes, a process is described for 
interference reduction by offset correcting a signal, where 
the offset may be modeled as including a step function and the 
signal is corrected by removing the undesired step function. 
For clarity of exposition, the description generally proceeds 
from an account of general elements and their high level 
relationship to a detailed account of illustrative roles, 
configurations, and components of the elements. 

Referring to Fig. 1, a generalized system 100 (e.g., a 
global system for mobile communications (GSM) , a time- 
division, multiple-access (TDMA) system, or a frequency- 
division, multiple access (FDMA) system) may be used to 
receive a transmitted signal and to correct the offset of the 
received signal, where the signal offset may be modeled as 
including a step function. Exemplary components of the system 
100 are described in greater detail below. 



The system 100 of Fig. 1 generally includes a transmitter 
110, a receiver 130 (e.g., a superheterodyne receiver, a dual- 
conversion superheterodyne receiver, or a direct conversion 
receiver) , and a channel 150 that models how the environment 
has changed the transmitted signal as perceived at the input 
port of the receiver. 

In general, the transmitter 110 and the receiver 130 may 
include any devices, systems, pieces of code, and/or 
combinations of these that may be used to transmit or receive, 
respectively, a waveform z(t) that generally may be 
represented as 

z(0 = Re {s{t)} cosicD^t) - Im {s{t)} sm(a)Qt) . (1.0) 
In equation (1.0), s(t) may denote a complex signal, coq = 2nfo 
may be an associated carrier frequency, and Re{s(t)} and 
Im{s(t)} denote respectively, real and imaginary parts of 
s(t) . 

A transmitter (e.g., transmitter 110) and/or a receiver 
(e.g., receiver 130) generally may include, for example, a 
mixer (e.g., mixer 135), a summer, a phase locked loop, a 
frequency synthesizer, a filter (e.g., filter 137), an 
oscillator (e.g., local oscillator 133), a frequency divider, 
a phase modulator, a down converter, an amplifier, a phase 
shifter, an analog-to-digital (A/D) converter or a digital-to- 
analog (D/A) converter (e.g., A/D converter 137), a 



microprocessor (MPU) , a digital signal processor (DSP) , a 
computer, or a signal processing circuit, whether linear or 
nonlinear, analog or digital, and/or any combination of these 
elements . 

More specifically, receiver 13 0 may include a down- 
converter for down- converting an input signal from radio 
frequency (RF) to baseband. The down- converter includes the 
local oscillator 133, the mixer 135, and the filter 137 (e.g., 
an infinite impulse response filter, a finite impulse response 
filter) . The receiver 13 0 also may include an A/D converter 
to generate a discrete signal from a continuous input and, for 
performing step parameter estimation and offset correction, 
any device, system, or piece of code suitable for that task, 
such as, for example, step parameter estimator and signal 
offset corrector 139. The step parameter estimator and signal 
offset corrector 13 9 may include, for example, a 
microprocessor control unit (MCU) , a digital signal processing 
(DSP) component, a computer, a piece of code, a signal 
processing circuit, whether linear or nonlinear, analog or 
digital, and/or any combination of these for use in performing 
the step parameter estimation and/or the offset correction, 
including step function removal. 

The transmitter 110 transmits a signal z(t) over the 
channel 150. The channel 150 may include any medium over 



which a signal may be communicated, such as, for example, an 
RF (radio frequency) portion of the electromagnetic spectrum, 
and or any other portion of the electromagnetic spectrum. 
Associated with the channel are an interference source 151 
that generates an interference signal J(t) and a noise source 
153 that generates a noise signal w{t) . The noise source 153 
and the interference source 151 add noise w{t) and 
interference J(t), respectively, to z(t) to form a signal r(t) 
received by the receiver. 

The noise w{t) may include, for example, additive white 
Gaussian noise that may have a zero or non-zero mean, while 
the interfering signal r(t) may have very different 
characteristics before and after an event that occurs within 
the burst. For example, if the interfering source is due to a 
different TDMA user who is transmitting at the same frequency 
as, but not time aligned with, the desired user, then the 
interference may appear as being turned on and off during the 
burst for multiple bursts. 

Fig. 2 illustrates one example of burst interference that 
may generate a signal offset that may be represented as 
including a step function at the input of the data bit 
estimator 140. The data bit estimator 140 may include, for 
example, a matched filter, and/or a decoder such as, for 
example, a convolution decoder, and may perform functions 
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including de-interleaving or decoding, and further may provide 
an estimate of data bits sent over the channel 150. 

Note that the signals z(t) and I(t), as shown in Fig. 2, 
are illustrative only and may not represent certain 
characteristics of actual physical signals. As shown, signal 
z(t) is transmitted in time slot n of a first TDMA waveform 
210, while the interference signal J(t) is transmitted in time 
slot m of a second TDMA waveform 220. Each TDMA waveform is 
associated with different TDMA channels (e.g., different GSM 
base stations with or without different hopping patterns) . 

The time slots for these first and second TDMA waveforms 
are not time aligned with each other (e.g., each time slot of 
the second TDMA waveform lags (or leads) the corresponding 
time slot of the first TDMA channel by the same time increment 
of t2-ti) . The interference signal J(t) also may have a power 
that is much greater than that of z{t) and a center frequency 
different than the center frequency coo of z(t), such as, for 
example, a center frequency that approximates a harmonic of 

COq. 

Because of the phase difference between the two TDMA 
waveforms, the transmission of signal J(t) at time t2 may 
appear as interference that is turned on and off and is 
included in signal z(t). Moreover, a TDMA channel structure, 
such as, for example, a TDMA time-slot assignment methodology. 



may ensure that z{t) and I(t) transmit in lockstep, causing 
z{t) to experience burst interference from I(t) beginning at 
the same relative point in each time-slot in which z(t) is 
transmitted (e.g., t2-ti from the beginning of each time slot). 

Referring again to Fig. 1, the receiver 13 0 receives from 
the channel a signal r(t) that includes z{t) + J(t) -h w{t) . A 
mixer 135 produces y(t) by mixing r{t) with the combination of 
a sinusoidal signal Q{t) generated by local oscillator 133 
with an attenuated version of r(t) that may leak into the 
local oscillator 133. The leakage of r(t) into the local 
oscillator 133 is represented by multiplying the received 
signal r(t) by an attenuation factor y to produce yr(t), and 
then summing yr(t) with the output Q{t) of an ideal local 
oscillator 134. Leakage of r(t) into the local oscillator 133 
causes y(t) to include the signal mix of r{t) [Q{t) + yr(t)] . 

Thereafter, y(t) passes through a low pass or band pass 
filter and/or an A/D (analog-to-digital) converter 137 (e.g., 
an integrator that performs the functions of A/D conversion 
and low pass filtering) to produce a discrete signal yn that 
may include an undesirable offset component. Thereafter, 
is processed further by step parameter estimator and offset 
corrector 139, which models the offset as a step function and 
estimates parameters descriptive of the step function. The 
offset of Yn is corrected by offset corrector 13 9 based on the 
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estimated step function parameters. 

Fig. 3 illustrates a receiver 130 that may be used to 
implement the system of Fig. 1, and in which a signal r{t) 
leaks into an ideal local oscillator 334. The signal r(t) may 
include a transmitted signal z{t) , a sum of interfering 
signals J(t), and additive white Gaussian noise signal w(t) . 
Signal z{t) may be represented as a real signal resulting from 
upconversion of a complex signal s(t) : 

zrt;=Re{s (t) }cos o)ot-Im{s (t) }sin coot. (1.1) 
P 10 Due to the leakage of r(t) into the local oscillator 333, the 

mixer 335 may not simply multiply r(t) by a sinusoid (e ,g . , AQe"^'""'^ , 
where Aq is a known value) . Instead r(t) is multiplied by the 
1*^ sinusoid and an attenuated version of the input, yr(t). The 

?! I 

resulting signal may be expressed as: 



0 
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y(t) = 4,[z(0 + 7(0 + w(t)y^''«' + r[\t) + ^\t) + yw\t) + 
2rz{t)I{t) + 2yz{t)w{i) + 2yw{t)I{t). 

Substituting equation 1.1 for z(t) produces 



(1.2) 



y{t) = 4,[Re{5(0}coscy/-Im{5(0}sinOo^]e-^'"°' 
+ 4>^(0e-^"°' + 4,w(0e-^"»' + ]d' (t) 

+ r[Re{s{t)}cos6)J - Im{5(0}sm co^tj + pv\t) (1.3) 

+ 2x[Re{5(0}cosfi>of - lm{s(t)}sm(o^t]r{t) 

+ 2/[Re{5(0}cos<yo? - lm{s(t)}smo)ot]Mt) + 2yw(t)I(t). 

If the attenuation term y is sufficiently small compared to the 
signal amplitude, then y(t) may be approximated as 
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y(t) » AQ\RQ{s{t)}coscoQt - Im{5(0}siniyQr]e" 



(1.4) 



in which the term yl^it) is retained because it is assumed that 
J(t) is of substantially greater power than z(t) . 

The signal y(t) then passes through a low pass filter 339, for 
example, to produce yiowit) , where yiowit) may be approximated as: 



in which the term AJ{t)e~^'^'^ is assumed to be substantially removed 
by the low pass filter and, therefore, has been dropped. In 
equation (1.5), the term w^^{t) represents a baseband portion of 



33 9. Assuming J(t) to be generally sinusoidal, yl'^{t) may include an 
offset component (e.g., a DC offset) and a bandpass component at 
twice the center frequency of r(t) . The low pass filter may 
substantially remove the bandpass component of //^(O while leaving 
essentially unaffected the offset component, represented in 
equation (1.5) as yll^ . When for example, J{t) is switched on or 
off, the offset component of yll^ may be modeled as a step function. 

An A/D converter 341 may be used to generate a discrete 
signal based on the signal yiowit) . Assuming that yll^ may 
be represented as a step function, the discrete signal yn may 
be represented as: 




(1.5) 



A^w{t)e 



that remains after passage through the low pass filter 
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5„(^) + cl + (c2-clK_^+>v„ 



(1.6) 



where Snid) is a discrete model of the baseband signal, 0 is a 
vector of unknown signal parameters (e.g., data bits, frequency 
offset) , and Wn is a discrete representation of zero-mean 
additive white Gaussian noise remaining after passing w{t) 
through the low pass filter 339 and the A/D converter 341. 
Also, referring now to Fig. 4, Un represents a unit step 
function that transitions from zero to one at n equals zero, 
such that cl + (c2 - d)M^_^ represents a step function with amplitude 
of cl before the step transition and amplitude of c2 after the 
step transition, where the step transition occurs at time n 
equals a. 

Referring again to Fig. 3, the signal is provided to the 
parameter estimator and offset corrector 139. The parameters 
cl, c2 and a of the step function are estimated, and the 
estimated parameters then are used to correct the offset of 
signal y„, to produce an output signal that may be represented 
as : 



The parameter estimator 13 9 may estimate the step function 
parameters based on, for example, gradient descent algorithms 
(e.g., the least -mean- square algorithm, Newton's method, the 



(1.7) 
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steepest descent method, and/or any combination of these methods) 
and/ oic the maximum- likelihood (ML) method. 

The ML method provides a general method of maximizing the 
likelihood of the joint probability density function of the 
values of the received signal vector {yo, . . . ,yN-i) given an 
intended signal vector (xo, . . . ,Xjj-i) . For the case when the 
observations are independent, a combined probability, or 
likelihood function, may be expressed as the product of the 
probability densities of the independent received signal vector 
samples, i.e., p = piyo)'-P(yN~i) , where it may be assumed that 
each probability density can be modeled as a Gaussian density. 
The likelihood function p is then maximized to find the optimal 
parameters using any suitable optimization technique (e.g., a 
non-linear optimization technique) , such as, for example, the 
Nelder-Mead method (a method based upon the simplex algorithm) , 
the steepest descent method, the LMS (least-mean-square) 
method, the Levenberg-Marquardt method (a least squares 
approach) , the Davidson- Fletcher- Powell method (a quasi -Newton 
based method) , or the Broyden-Fletcher-Goldf arb-Shannon method 
(a quasi-Newton based method) , and/or any combination of one or 
more of these or other optimization methods. 

More specifically, a ML estimate of the step function 
parameters cl, c2, and a can be obtained from the samples as 

described above. For example, we may take the baseband signal 
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model s^(6) and noise model w(n) to have a zero mean, since 
their means can be incorporated into the step function 

parameters. Using -^s^(d) + cl + (c2-cl)u^_^ as an expression of the 

mean of the individual values of the received signal vector 
produces the following ML likelihood function of the complex 
observation: 

„=o V^cr 

which may be simplified to 



P = 



( 1 \^ -Vi;L~^«(^)+^^i+(^2-ci)«^| 

^ ^ '^'^ ' . (1.9) 



To maximize the value of p, it is sufficient to minimize 
the value of 



«=0 



yn-^sM-c\-ic2-c\)u\' , (1.10) 



2 

which is a nonlinear least squares optimization problem. 
Specifically, the unknown parameters in equation (1.6) can be 
determined by solving 



N-\ 

min /* = y 

n—0 



y,-^sM-c\-{c2-c\)u,^^ ' . (1.11) 



2 

To determine the solution of (1.11), it is useful to 

partition equation (1.11) over a first interval before the 

transition of the square wave and a second interval after the 

transition of the square wave. That is, equation (1.11) 
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becomes : 



n=0 



N-\ 



+ 



(1.12) 



Equation (1.12) may be minimized over {9, cl, c2, a) 
jointly using any of the previously mentioned optimization 
methods. However, since cl and c2 are in separate terms of the 
objective function, their estimates also may be solved for 
separately. For example, the estimate for cl may be obtained 
analytically by differentiating the portion of equation (1.12) 
that corresponds to the first interval with respect to cl, 
setting the result equal to zero, and solving for cl . The 
estimate of c2 may be solved by operating upon the portion of 
equation (1.12) that corresponds to the second interval in like 
fashion. 

The estimates of cl and c2 also may be obtained 
qualitatively. For example, the estimate for cl may be 
expressed as a mean of an error between the observation and 
the signal prediction Snid) before the square wave transitions; 
similarly, the estimate for c2 may be expressed as a mean of an 
error between the observation and the signal prediction SniO) 
after the square wave transitions. Hence, the estimates cl of 



cl and c2 of c2 may be expressed as 



1 a-l 
^ n=0 



(1.13) 
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and 



1 N-l 

c2=-L-y 



(1.14) 



Equations (1.13) and (1.14) then may be substituted back 
into the objective function of equation (1.12), resulting in 
the following expression of the objective function: 



a-l 



J 

N-\ 



1 a-l A 1 a-l A 



m=(i 
N-a 



N-a 



(1.15) 



Now, equation (1,15) is a function of the obseirvation, yn, 
the unknown signal parameters, 0 , and the location of the step 
function, a. All of these parameters may be jointly estimated, 
for example, using non-linear optimization techniques as 
described above. 

Nevertheless, it also may be possible to estimate only the 
unknown parameters cl, c2 and a based on expanding and 
rearranging the terms of equation (1.15) to give 



n=0 



-a 



« «=o ^ 



1 " 

m=0 



+ 2aRe 
+ 2(iV-a)Re 



-(N-a) 

1 N-a 

l^y,. 



1 ^^Ai 



(1.16) 



1 ^^"4, 
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X 



where 



g(a) = a 



1 ^ 



1 



(1-17) 



The first term in equation (1.16) is an expression of mean 
square error between the observation yn and the signal 

prediction -^s^Gi) , while the second term is explicitly provided 

fc'^ by equation (1.17) . 

All of the other terms of equation (1.16) involve averages 



|1 of the signal prediction — 5^(^)and may be approximated as zero 

^ 10 if the expectation of s^{6) is approximately equal to zero, for 

m 

both before and after the transition of the step function. When 

%i 

23 the expectation of s^{6) may not be approximated as zero, the 



parameters may be estimated using a method that retains these 
terms. For example, the parameters may be estimated by starting 
15 with a seed value of a that may be used to determine an 
estimate of 9 , which, in turn, may be used to produce an 
estimate of a. The method may be iterative and may continue to 
alternate between estimating a and 9 until convergence to a 
desired degree of precision is achieved. 

Nevertheless, for many signals, such as, for example, a GSM 
signal for which the expected value of the underlying binairy 
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m 



%4 



15 



20 



data stream is zero or approximately zero, it is reasonable to 
assume that the average of the signal prediction Sn{6) is equal 
or approximately equal to zero when taken over a sufficiently 
large interval. For example, the signal prediction may be 
expressed as : 



(1.18) 



A:=0 



where djc is an original binary data sequence with an expectation 
of zero, and hn is the combined action of the transmit filter, 
the channel filter, and the receive filter. Here the vector of 
unknown parameters, 0 , can be taken as the complete data 
sequence djc for all k and the complete filter hn for all n. 
Because the expectation of the binary sequence is zero and the 
binary sequence is independent of the combined filter, then the 
expectation of the signal in equation (1.18) is zero. That is, 

^K} = 0 implies E{s^(9)} = 0 ^ (1.19) 
Hence, approximating as zero the expectation of Sn{Q) , the 
objective function of (1.16) becomes: 



n=0 ^ 



(1.20) 



and equation (1.20) may be minimized by selecting an cx that 
maximizes g{a) . That is. 



a » argmaxg(a) = argmaxa 



1 V 

-Z^yn 



+ (N-a) 



1 

N-a 



i\ —u. 



(1.21) 
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Fig. 5 illustrates a method 139 for optimizing equation (1,21) 
that may be used to implement the system of Fig. 1. A sum of the 
received data is computed and stored (step 510) , where the sum may 
be expressed as: 



Y,-Y,yn- (1.22) 

B=0 



Next, temporary parameters Yps (a partial sum of the data) 
and graax are set initially to zero, and temporary parameter ^rest 
is set initially equal to one (step 520) • 

Using the parameters of step 520, estimates gmaxr ol , and Y^^ 

may be computed iteratively over increasing values of a-r&st while 
a^est is less than or equal to N-1, the number of data samples 
(steps 530) . For example, as shown in Fig, 5, estimating gmax/ 

d , and Y^^ may include adding the current data Y^^ to the 

partial sum of the data Yps to generate an updated partial sum Yps 
(step 533) . An updated value for the object function g then may 
be determined as : 

(step 535) . The updated value of g then may be compared to the 
stored value of g^aK (step 537) , and if updated g is greater than 
gmaxi then gr^nax tnay be set equal to updated gr as a best current 

guess of the maximum of g, a may be set equal to arest/ and Y^^ 

may be set equal to Yps (step 539) • After updating the values of 
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gmaxi OL f and 7^^ {step 539), a^est tnay be incremented (step 541) 
and, if arest is less than or equal to N-1 (step 531) , then steps 
530 may be repeated. 

The estimation of the parameters accomplished in steps 53 0 
also may be performed, for example, by decrementing a^est from a 
high value to a low value, or by performing a random selection 
of aTest- Under any of the approaches mentioned, parameters may 
or may not be estimated for each value of ofrest- 

Following completion of the iterative process of steps 530, 

the estimated values of gj^ax, a, and Y^^ may be used to correct 

the offset of the data y„ (step 550) . For example, using a as 
the estimate of the transition point of the step function, the 
estimate c\ may be expressed using the calculated values as 
1 - 

^1 = -^;./ (1.24) 
a 

while c2 may be expressed as 

^2 = ^(7,-i;J. (1.25) 

Optionally, where c\ and c2 as estimated above are equal or 
approximately equal (indicating that a step function may not be 
present) , then both may be re-estimated as 

cl = c2 = l(yj. (1.26) 

Thereafter, to correct the offset of the received data y^, the 
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estimated parameters may be used to subtract the step function 
from each data point as follows 



Following the correction of the offset, further estimation 
methods may be applied to the residual data (y^ minus the step 
function) in order to estimate the remaining unknown signal 
parameters 6. 

Other implementations are within the scope of the 
following claims . 




Q<n<a 



(1.2 
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